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We investigate the phase diagram of an extended Kondo lattice model for doped manganese oxides 
in the presence of strong but finite Hund's coupling and on-site Coulomb interaction. By means of 
the Schwinger-boson mean-field approach, it is found that, besides magnetic ordering, there will be 
non-uniform charge distributions, such as charge ordering and phase separation, if the interaction 
between electrons prevails over the hybridization. Which of the charge ordering and phase separation 
appears is determined by a competition between effective repulsive and attractive interactions due 
to virtual processes of electron hopping. Calculated results show that strong electron correlations 
caused by the on-site Coulomb interaction as well as the finite Hund's coupling play an important 
role in the magnetic ordering and charge distribution at low temperatures. 



I. INTRODUCTION 

The doped manganese oxides with perovskite struc- 
ture, Ri-^A^jMnOa (R=La, Pr, Nd; and A=Sr, Ca, 
Ba, Pb), have recently attracted much attention duejntpi 
the observation of colossal magnetoresistance (CMR) 
Theories based on the double exchange (DE) model, 
in which the exchange of electrons between neighboring 
manganese ions with a strong Hund's coupling drives 
spins of on-site electrons to align parallelly, have been 
developed for a long time and can qualitatively elucidate 
the relation of transport and magnetism in the doping 
range of 0.2 < x < 0.45. However, recent systematic ex- 
perimental studies revealed rich phase diagrams, which 
are difficult to be understood by the DE model alone. 
For instance, the system is actually insulating in the half- 
doping case {x = 0.5) at low temperature; but a metallic 
ferromagnetic (EM) state would be expected according 
to the DE model, for the DE hopping reaches its maxi- 
mum and the effective EM interaction becomes strongest 
at this doping. Eurthermore, for x = 0.5, a charge order- 
ing characterized by an alternating Mn'^+ and Mn^+-ion 
arrangement in the real space was observed to be 
perimposed on the antiferromagnetic (AE) ordering 
This charge ordering is sensitive to an applied magnetic 
field, and even melts under a moderate magnetic field. 
In the meantime, thejjcsistance may decrease by several 
orders in magnitude,Llil3 implying that there is a close 
relation between the charge ordering and the AE spin 
background. 

Many efforts have been made to understand the phase 
diagram_nL^ doped manganese oxide based on various 
models.EirLj To explore the origin of the unique mag- 
netotransport, many further theories on the basis of 
DE model have been proposed, such as Jahn-iTellfir dis- 



placement and elep 
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■ron-phonon interactions 
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nal and off-diaffiaaaJ disorder,E3L3 and the phase sepa- 
ration scenario.EjS Importance of various interactions 
on the physical properties of the manganites is currently 
one of the lively-debating subjects. A comprehensive un- 
derstanding of the magnetic and charge ordering states 
as well as their relations to the transport properties are 
highly desirable. 

In this paper, we investigate effects induced by strong 
but finite Hund's coupling Jh and on-site Coulomb in- 
teraction U . In the manganites, three t2g electrons are 
almost localized and form an S" = 3/2 spin state accord- 
ing to the Hund's rule, while electrons in Cg orbit evolve 
a conduction band. In the conventional DE model Jh is 
usually regarded to be infinite and U is neglected, so that 
there exist only single-occupied state with spin 5+1/2 
and empty state with spin S. In this limit, since neither 
S" — 1/2 state nor double occupancy of Cg electrons on the 
same site is allowed, the effects induced by two kinds of 
virtual processes, that an electron hops to a empty site 
with its spin antiparallel to the core spin on the site and 
that an electron hops to a singly-occupied site, are com- 
pletely eliminated. In this work we will show that the 
effects due to the two kinds of virtual processes are im- 
portant to account for the observed phase diagrams. On 
one hand, the two kinds of virtual processes can produce 
an AE superexchange coupling between the neighboring 
localized spins, with coupfing strength being i^/(2J/f) 
and / {Jh + U), respectively, where t is the hopping in- 
tegral of the eg electron. The AE coupling induced by 
the virtual processes is usually much stronger than the 
direct AE superexchange coupling between neighboring 
localized spins. On the other hand, the virtual process 
of producing a single-occupied state with spin — 1/2 
can lead to a repulsive interaction between conduction 
electrons, while the virtual process of producing a dou- 
ble occupancy can result in an attractive one, as will be 
shown in Sec. HA. The two types of interactions between 
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electrons can produce non- uniform charge distributions, 
such as charge ordering and electronic phase separation, 
provided that the interactions prevail over the hybridiza- 
tion effects of electrons. Which of charge ordering and 
phase separation appears is determined by a competition 
at different doping between two types of interactions, for 
the repulsive interaction is favorable to the charge order- 
ing while the attractive interaction may cause the phase 
separation. 

In Sec. II, starting from an extended Kondo lattice 
model, we have an effective projected Hamiltonian in the 
case of strong but finite Jh and U. The repulsive and 
attractive interactions, which are associated with mag- 
netic ordering and non-uniform charge distributions, are 
obtained. A Schwinger-boson mean-field theory is devel- 
oped to establish the phase diagram at low temperatures. 
In Sec. Ill we focus attention on the case of a; = 1/2. 
The possibility that the FM, AF or canted ferromagnetic 
(CF) order appears, as well as that the Wigner lattice 
or phase separation appears, is discussed. In Sec. IV 
numerical results for phase digram are presented. In the 
large case, the hybridization effect is dominant and 
there exists metallic ferromagnetism, which accords with 
the DE model. As the repulsive interaction due to fi- 
nite Jh effects is relatively strong, the charge ordering 
will be formed at half doping in the AF background. On 
the other hand, the phase separation may arise when the 
attractive interaction due to finite U becomes dominant. 



II. GENERAL FORMALISM 

A. Effective projected Hamiltonian 

The electronic model Hamiltonian for doped man- 
ganese oxides we considered presently is given by 

H = -tY^ cj^Cj^ + U 

{ij},cr i 

-Jh ^ cL^J • CTcrcr'Cio-' + J AF ^ • Sj . (1) 

iaa' (ij) 

Here Cia {c]^) is the annihilation (creation) operator for 
conduction electrons at site i with spin a, ni^ = cl^Cia is 
the particle number operator, Si is the total spin operator 
of the localized electrons at site i with S = 3/2, and cr is 
the Pauli matrix. In the manganitcs the Hund's coupling 
constant Jh and the on-site Coulomb interaction U are 
much greater than the hopping integral t as well as the 
direct AF superexchange coupling Jaf- On a given site, 
an itinerant electron constrained by the strong Hund's 
coupling has its spin parallel to the core spin, forming a 
spin 5* -I- 1/2 state. The singly-occupied state for the itin- 
erant electron with spin antiparallel to the core spin and 
the doubly-occupied state are almost forbidden, making 
it appropriate to utilize the projective perturbation tech- 
nique to investigate Hamiltonian (1). The effects of finite 



Jh and U can be regarded as a perturbation correction 
to the large Jh and U limit where there are only the 
empty and single occupancy with S* + 1/2 state. Up to 
the second-order perturbaticp, an effective Hamiltonian 
for Eq. (1) can be derived aj£3 

^^ff = H ^^Cj'^ + Jaf Si ■ Sj 

{ij),<^ (ij) 

+j,Yi^,-s,-smhPi 

(y) 

+J2j2i^^-S,-S')PlPl . (2) 
(y) 

Here Si is the spin operator with S ^ S + 1/2. Pui 
and Pjt, are the projection operators for empty state with 
spin Si and single occupancy with spin Si, respectively. 
S, = S,P,h + SS,pI/S. = P,c,,P, and 4 = P^clP, 
are projected electron operators, where Pi = Pih + Pj^, 
projects onto the space of non-doubly-occupied site. The 
last two terms in Eq. (2) are the second-order pertur- 
bation corrections where Ji = /{2JhS^) and J2 = 
/[{JhS + U)S'^], respectively, stemming from differ- 
ent virtual processes (a) and (b). In virtual process (a), 
an electron first hops from a site to one of the nearest- 
neighbor empty sites to form a spin S — 1/2 state and 
then hops backward; while in virtual process (b), an elec- 
tron first hops to a single-occupied site, where there has 
been an electron with opposite spin, and then hops back- 
ward. Owing to Ji > and J2 > 0, both virtual pro- 
cesses favor the AF arrangement of the core spins and 
enhance greatly the AF coupling between the neighbor- 
ing spins. Whether non-uniform charge distributions can 
appear is dertermined by a competition between elec- 
tronic hybridization and interactions of electrons with 
one another. The hybridization, the overlap of electron 
wavefunctions centred on different sites, is a quantum- 
mechanical effect that allows electrons to hop from one 
atom to another, thus tending to spread the electronic 
density uniformly through the system. In contrast, the 
interactions of electrons in the sjstem tend to promote 
non-uniform charge distributions.tfl The addition of the 
last two terms in Eq. (2), arising from the finite Jh and U 
effects, enhances greatly the interaction side in the com- 
petition and thus favors the occurance of the non-uniform 
charge distributions. On the other hand, since the val- 
ues of Ji{Si ■ Sj — SS) and J2(Si • Sj — S^) are always 
non-positive, it is clear that the two terms represent the 
repulsive and attractive interactions between conduction 
electrons on neighboring sites, respectively. The compe- 
tition between them will lead to different non-uniform 
charge distributions, the charge-ordered state and the 
phase separation. 
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B. Schwinger boson representation 

A projected electron operator may be regarded as 
a combination of the operator / for a spinless charge 
fermion and that for a neutral boson with spin S = 

5 + 1/2, i.e., Cia rnfibicr/ where bia is the Schwinger 
boson operator Bj In the Schwinger boson representa- 
tion, the spin (S'i = Si or Si) can be expressed as 
= hT^aa'^L'^'^'^'^'^'^'' with constraint J2a4abta = 
2S' {S' = S oi S). This constraint can also be written 

as J2cr ^icr^icr = 2S + Ui, whcrc rii = fj fi is the particle 
number operator of fermions. The projection operators 
can be replaced by Pit = 1 — ni and P^, — rii. Then, 
by using the identity S', • S^ = -\A\jA,j + S',S'j with 
■^ij ~ hibi]bji — biibji), the effective Hamiltonian is re- 
duced to 

H^JAFS'-iJ2fIfjP^J 

-^AljAij[JAF + Jii^ - ni)nj + J2ninj] , (3) 

where F,^ = EAI^j<t, i =t/{2S), Ji = Ji-Jaf/S and 
J2 = J2 — Jaf/S + Jaf/^S"^- The difference between J„ 
and J„ (n =1 and 2) is caused by the finite S effect, and 
would disappear in the large spin limit. 

C. Mean-field approximation 

We now make a Hartree-Fock mean-field approxima- 
tion by decoupling various terms in Eq.(3). As an exam- 
ple, we have 

AljAijUiUj AIj < Aij X >< rij > 

+ < Alj > A,j < Ui >< rij > 

+ < Alj >< A,j > n,< rij > 

+ < AIj >< Aij >< rii > rij 

-3 < A\j X Aij >< rii >< rij > . 

For the fermions, we have < rii 1 — a; in a uniform 
density state. In the presence of the charge ordering at 
X = 1/2, the lattice can be divided into two sublattices 
X — A and B, and < rii >— Nx = 1 — x + aS where 
Q = 1 or — 1 for i E A or B, and 5 (0 < (5 < 1 — x) is 
the charge ordering parameter. The uniform state.. can be 
regarded as a special case of S — 0. As in Ref.t3£3, the 
FM and AF order parameters can be respectively written 
as 

<F..,> = F , 

< A.,, > ^ i-ipA , 

where ri depends on the position of site i, being an odd 
number in one sublattice and an even number in the other 



sublattice. The average value of < fjfj > is assumed to 
be a constant K. Under the Hartree-Fock mean-field 
approximation, Hamiltonian (3) is reduced to 

Hmf — Eq + Hbosc + Flpermi (4) 

with 

Eo ^ NZiKF + NZJafS^ - NA2S 
+ ^NZA^Jaf + 2(1 - x).h 
+3[il - x f - S%h - Ji)} 

ia- (ij)a 

Hpernii = -tF ^ j} f j - (A + ^) ^ ni 

(y) i 

-]^A^zY}%h ~ Ji)Nx + Ji]n, ■ 

i 

Here J= JAF + {J2-Ji)[{'f-xf-6'^] + Ji{l-x) is the 
effective AF coupling constant, N is the number of lattice 
sites, Z is the coordinate number of the lattice, and N-^ = 
Nb {Na) for i e ^ {B). A =< Ai > with A^ being 
the Lagrange multiplier introduced at site i to enforce 
the local constraint b\^bi„ = 25-1- riif^ and jj. is the 
chemical potential of fermions which is determined by the 
total conduction electron number. After diagonalizing 
Huose and Hp emu, we have 

Huose = ^k/^L/^k. + ^(Pk - A) , 
kcr k 

and 

H Fermi = ^(ek - /i)V'I'/'k , 
k 

with boson and fermion dispersions as 

= ^A2-(ZAJ7k)2 - Zt7^7k , 

and 

Ek = -S5n(7k)Z^A4(J2 - Ji)2^2 + I2p2^2 ^ 

respectively. Here 7k = Z^^ J^r/ e'*''^ with rj the vec- 
tor of the nearest neighbors of each site, and the shifted 
chemical potential is given by 

/i = /i + A + ^ZA^[Ji + 2(J2 - Ji)(l - x)] . 

From the fermion spectrum, it is found that the energy 
band of fermions is divided into two separated branches, 
between which there is an energy gap provided (5 7^ and 
J2 7^ Ji- Thus, in the half-doping case {x = 1/2) at zero 
temperature, the lower branch (7k > 0) is fully occupied 
by fermions while the upper one (7k < 0) is empty. 
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III. CHARGE ORDERING AND PHASE 
SEPARATION 



with i = 1,2. Then the set of mean- field equations (5a)- 
(5d) can be rewritten as 



Now we come to discuss the magnetic ordering and 
non-uniform charge distributions at a; = 0.5. First let us 
write down the mean-field equations. When the magnetic 
ordering arises, the Schwinger bosons should condensate 
to the lowest energy state. From £'k=o = 0, we have 



A = zV¥io+A^ . 

Using the spectra of quasi- fermions and bosons, we ob- 
tain a set of self-consistent equations at zero tempera- 
ture: 



F = 



2^ ~ ^(Ji - J2)^A^^ + t^F^^l 
iK{2S + 2-x) 



(5a) 



tK 



A 



^i-^K^ + A^.P{l--fl) 
AJ(2S + 2-x) 
Vi^K^+A^P 



(5b) 



N 



k .^i^K^ + A^P{l--tl) 

O N 2^ 



k J(Jl - J2)M4(52+t2^2^2 



(5c) 
(5d) 



From Eq. (5d), it follows that the charge ordering may 
appear {6 ^ 0) only when Ji > J2, otherwise there exists 
only a trivial solution oi S = 0. Physically, the effective 
interaction between electrons is repulsive for Ji — J2 > 0, 
which is responsible for the formation of the charge or- 
dering. There are two types of possible ordering states: 
the pure magnetic ordering phase without charge order 
and the charge ordering phase in which the charge and 
magnetic orders coexist. The former includes FM, AF, 
and CF phases; while the latter is the Wigner lattice 
{6 = 1/2, K = F = 0, and A = A^ax) or a combina- 
tion of the charge ordering (0 < 5 < 1/2) and magnetic 
ordering {A < A^ax, K and F ^ 0). Since the system 
under consideration is a three-dimensional (3D) simple- 
cubic lattice, for convenience, wc introduce two param- 
eters m = PfV[(Ji - J2fA'^S^] {0 < m < 0°) and 
= -A^P/{PK^ + A^.P) (-1 < 7^2 < 0) and define 



the following two 3D integrals: 



K 



771D1 



F = V1+^(25 + 2-X-C2) , 
A = V^C^S + 2-X-C2 + D2) 



(6a) 

(6b) 
(6c) 

(6d) 



For the Wigner lattice {6 =1/2), from Eq. (|^), it fol- 
lows Ci = 1 and ?/i = 0, so that K — Q, rj2 = — 1, 
and _F = 0. In this case, A reaches its maximum 
Amax ~ 2S + 0.597, indicating a fully AF insulator in 
which all the electrons are localiezed in one sublattice. 
Another limiting case is that there is no charge order- 
ing {5 ~ 0), where there are two sets of solutions: (i) 
K = F = and A = Amax, which is an AF state; and 
(ii) rji — -foo, K ^ Q, and F ^ Q, which is a CF state 
for < ^ < Amax or a FM state for A = Q. There is 
the lowest energy in the ground state. This ground state 
energy per site in the AF state [case (i)] is 



J'^maxl'^ 



and that for case (ii) is given by 



{2S + 2-X-C2- mD2)i 
%^J1 + ri2 



Pri2 



72J(l + r;2) 



(7) 



(8) 



which is a function of 772. For a CF state, the magnitude 
of •q2 is determined by minimizing Eg, i.e., dEg/dr]2 = 0. 
7^2 = corresponding to a FM state. It is found that a 
phase transition occurs from a CF state to a FM state 
when J decreases to t/{12S + 9). 

The uniform density phases discussed above may be 
unstable toward the phase separation, which will occur 
if d^/dx < 0. In our model, when the virtual process 
(b) is dominant over process (a), i.e., J2 > Ji, the effec- 
tive interaction between the Cg electrons on neighboring 
sites will be attractive. If such an attractive interaction 
is strong enough, the Cg electrons will tend to accumu- 
late together to lower the total energy of the system. In 
this case, the electronic density is no longer uniform and 
the phase separation occurs between electron-rich and 
electron-poor regions. The chemical-potential judgment 
of phase separation is d^/dx < 0, equivalent to the en- 
ergy judgment d'^Eg/dx^ < due to /i = dEg/dx. The 
phase separation will appear when d^Egjdx^ < 0, where 
both r]2 and J in Eqs. (0) and (||) are functions of x. 



IV. PHASE DIAGRAMS AND DISCUSSIONS 

Mean- field equations (6a)- (6d) have been numerically 
solved and the phase diagrams at a; = 0.5 are plotted 
in Figs. 1 and 2. To make the results comparable with 
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each other for diflFerent S, we have used reduced cou- 
pling constants jaf = JafS"^ and jh = JhS. Fig. 1 
is the phase diagram of jaf/t and t/jh in the large U 
limit with S = 3/2 and S oo. In this case, J2 — 0, 
and t/jh indicates the ratio of the repulsive interaction 
due to virtue process (a) to the kinetic energy of the Cg 
electron. It is found that the system is in the FM state 
for small jaf/t and t/jh, where either direct or indirect 
AF superexchange coupling is weak compared with the 
FM coupling due to the DE mechanism. With increas- 
ing jaf/t and t/jh, the AF superexchange coupling be- 
tween the neighboring localized spins is enhanced. As 
J > t/{12S + 9), the CF state appears; meanwhile, a re- 
pulsive interaction between fermions is increased due to 
the increase of t/jh- At this stage, there is a CF order- 
ing but the charge ordering has not yet arisen, which is 
labeled as the CFl state in Figs. 1 and 2. When t/jh is 
increased beyond a threshold so that the repulsive poten- 
tial between fermions dominates over the kinetic energy, 
the fermions begin to tend toward one sublattice and 
thus charge ordering emerges, which is called as the CF2 
state. Finally, with further increase of t/jh, S increases 
gradually. At ^ = 1/2, all the Cg electrons are confined in 
one sublattice in the AF background, forming a Wigner 
lattice. It is found that the phase diagram for S* = 3/2 is 
very similar to that for S —^ 00. The difference between 
them can be seen in the large jaf/t case, for the effective 
repulsive interaction is proportional to the factor Ji — J2 
or t/jh— jaf / {2tS^S). For finite S, the increase of jaf will 
decrease this factor and so be unfavorable to the charge 
ordering. As a result, the quantum fluctuation of finite 
S leads to enlargements of the CFl and CF2 regions. 

For manganites, it is roughly estimated that t £i 
0.15eV, jh « 0.75eV, U w lOeV, and jaf w SmevjB 
so that jaf/t « 0.05 and t/jh « 0.2. From Fig. 1, it 
follows that the repulsive interaction from virtual pro- 
cess (a) alone enables the Cg electrons to form the charge 
ordering. This indicates that this process plays an impor- 
tant role in determining the collective behavior of elec- 
trons. In actual doped manganites, there may also be 
other effects that favor charge ordered phase, which are 
not considered in the present paper, sjjfhjas. the direct 
nearest- neighbor Coulomb interactions .EJ'EZrta It is worth 
pointing out here whether the charge ordering occurs de- 
pends not only the ratio of t/jh, but also the magnetic 
ordering of the system. In doped manganites, the ampli- 
tude for an electron to hop from one site to another is 
determined by the relative orientation between the core 
spins at the two sites, being greatest when the core spins 
are parallel and least when they are antiparallel. As a 
result, a FM state has the greatest hybridization and so 
tends to spread the electronic density uniformly through 
the system; while an AF state has the least hybridiza- 
tion and favors the charge ordering. This can account 
for the sensitivity of charge ordering with respect to an 
applied-magnetic field that tends always to align the core 
spinsB'Ej Besides, there is a difference between two types 



of repulsive interactions mentioned above in the response 
to an applied magnetic field: with a decrease in the AF 
correlation, the superexchange electron-electron interac- 
tion is reduced, while the direct Coulomb interaction re- 
mains unchanged. Therefore, charge ordering caused by 
virtual process (a) is more easily affected by the magnetic 
field than that caused by the direct Coulomb interaction. 





WL 


(a) 


\ \ CF2 






FM \ CFl 


\ \ '^^^ 


WL 


(b) 


FM \ 









FIG. 1. Phase diagram of the extended Kondo lattice 
model with [/ = 00 for (a) S = 1.5 and (b) S ^ 00. FM, 
CFl, and CF2 denote regimes of metallic ferromagnet, canted 
ferromagnet without and with charge ordering, respectively. 
WL stands for the Wigner lattice. 

In Fig. 2, we present the phase diagram of jh/ {jh + U) 
and t/jh- In the finite U case, jh/{jh + U) describes the 
relative strength of the effective interaction due to vir- 
tual process (b) to that due to process (a). For small 
t/jh, there is only a FM state, where the Hund's cou- 
pling is very strong, both virtual processes (a) and (b) 
are suppressed so that the hybridization effect and the 
DE ferromagnetism dominates over the system. With 
increasing t/jh, both Ji and J2 become large and so the 
AF superexchange coupling due to virtue processes (a) 
and (b) is enhanced, leading to the CFl state. As t/jh 
is increased beyond a threshold, the charge ordering or 
phase separation may occur, depending on the compe- 
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tition between effective repulsive and attractive interac- 
tions. For jh < U, Ji > J2, the net repulsive interaction 
favors the charge ordering. The opposite case is Ji < J2, 
where there is is a net attractive interaction and so the 
phase separation may occur. In the middle region near 
jh/{jh + U) = 0.5, a cancellation of Ji and J2 leads to a 
very small net interaction so that the non- uniform charge 
phase can not be formed. Thus, the phase diagram shown 
in Fig. 2 is determined by two types of competitions: one 
competition between hybridization and interaction and 
the other between the repulsive and attractive interac- 
tions, respectively, arising from virtual processes (a) and 
(b). 




i,/(j,+u) 



FIG. 2. Phase diagram of t/jh and jh/U. The AF cou- 
pling constant is taken to be jaf/t = 0.03. PS stands for the 
phase separation. 

At this stage, we wish to point out that in realistic 
manganites, the charge ordering is likely to be accompEu 
nied by the orbital ordering and the lattice distortion,^ 
which is not considered in the present model, for the main 
purpose of this work is to examine the effects of the fi- 
nite Jh and U. More virtual processes will be involved if 
the orbital degeneracy is taken into account, since there 
are more mediate states. Meanwhile the lattice distor- 
tion will affect the virtual processes due to its effect on 
the electron hopping. These additional virtual processes 
may be important in the realistic manganites. 

Finally, the present method is not confined to the half- 
doping case. Fig. 3 is the phase diagram for arbitrary 
doping, in which, for simplicity, U is set to zero to max- 
imize the effect of process (b). In this case, it is found 
that the phase separation occurs at either high or low 
concentrations of the Cg electrons. The phase diagram 
obtained hese,is very similar to the result of Monte Carlo 
simulations .Eil Similar-xesiilts were also obtained analyti- 
cally by other authorsESES. As U is increased, process (b) 
will be suppressed, since the energy paid necessarily for 
the double occupancy, U + JhS, will become high. The 
phase separation disappears when U is large enough. In 



doped manganites, the on-site Coulomb interaction U is 
usually much stronger than the Hund's coupling JhS. 
However, if the orbital degeneracy at each site is further 
taken into account, the energy paid for the doubly occu- 
pancy might be much lower than U + JhS, for two Cg 
electrons at the same site can occupy two different orbits 
and so their spins can be parallel to each other as well 
as the core spin. In this sense, the virtual process of the 
double occupancy may revive. It is worth mentioning 
that the coincidence of our result at C/ = with the nu- 
merical simulation indicates the reliability of our method 
in this system. 




0.0 0.2 0.4 O.B 0. 



1-x 



FIG. 3. Phase diagram for jaf/t = 0.01 and U = 0. 

V. CONCLUSION 

In summary, we have studied an extended Kondo lat- 
tice model in the presence of strong but finite Hund's 
coupling and on-site Coulomb interaction. By means of 
the Schwinger-boson representation and a mean-field ap- 
proximation, we show that the effects of the finite Jh and 
U favor antiferromagnetism. In the half-doping case, it is 
found that the charge ordering may be superimposed on 
the magnetic ordering by the repulsive interaction be- 
tween electrons due to the virtual process of electron 
hopping. If the on-site Coulomb interaction is weak, the 
phase separation may occur in either high- or low-doping 
case. The calculated results show that the finite Jh and 
U effects play an important role in forming the mag- 
netic ordering and non-uniform charge distributions in 
the doped manganese oxides. 

ACKNOWLEDGMENTS 

This work was supported by a CRGC grant at the Uni- 
versity of Hong Kong. One of us (D.Y.X) acknowledges 
support from the National Natural Science Foundation 
of China. 



6 



^ R.von Helmolt, et al., Phys. Rev. Lett. 71, 2331 (1993). 
^S. Jin, et al., Science 264, 413 (1994); M. McCormack et 

al., Appl. Phys. Lett. 64, 3045 (1994). 
^Y. Tokura et al., J. Phys. Soc. Jpn. 63, 3931 (1994); A. 

Asamitsu et al.. Nature (London) 373, 407 (1995). 

* J. Fontcubcrta et al., Phys. Rev. Lett. 76, 1122 (1996). 
^ C. Zener, Phys. Rev. 82, 403 (1951). 

^ P. W. Anderson and H. Hasegawa, Phys. Rev. 100, 675 
(1955). 

^P. G. de Gennes, Phys. Rev. 118, 141 (1960). 

* J. B. Goodcnough, Phys. Rev. 100, 564 (1955); K. Kubo 
and N. Ohata, ,J. Phys. Soc. Jpn. 33, 21 (1972). 

"H. Kuwahara ct al.. Science 270, 961 (1995); T. Vogt et 

al., Phys. Rev. B 54, 15303 (1996). 
" Y. Tomioka et al., Phys. Rev. Lett. 74, 5108 (1995). 
" J. Inouc and S. Maokawa, Phys. Rev. Lett. 74, 3407 (1995); 

M. Yamanaka, W. Koshibae and S. Maekawa, Phys. Rev. 

Lett. 81, 5604 (1998). 

R. Shiina, T. Nishitani and H. Shiba, J. Phys. Soc. Jpn. 
66, 3159 (1997). 

S. Ishihara, M. Yamanaka, and N. Nagaosa, Phys. Rev. 
B56, 686 (1997). 
^'^ J. van den Brink and S. Khomskii, Phys. Rev. Lett. 82, 

1016 (1999). 

D. P. Arovas and F. Guinea, Phys. Rev. B 58, 9150 (1998). 
A. J. Millis, P. B. Litlcwood, and B. I. Shraiman, Phys. 
Rev. Lett. 74, 5144 (1995); A. J. Millis, B. I. Shraiman, 
and R. Mueller, ibid. 77, 175 (1996). 

A. S. Aloxandrov and A. M. Bratkovsky, Phys. Rev. Lett. 

82, 141 (1999). 

T. T. M. Palstra ct al., Phys. Rev. B 56, 5104 (1997). 
" L. Sheng et al., Phys. Rev. Lett. 79, 1710 (1997); E. Miiller- 

Hartmann and E. Daggoto, Phys. Rev. B 54, R6819 (1996); 

C. M. Varma, Phys. Rev. B 54, 7328 (1996), D. Y. Xing, 

S. J. Xiong and D. Feng, Phys. Rev. B 58, 14139 (1998). 
2" F. Zhong, J. Dong and Z. D. Wang, Phys. Rev. B 58, 15310 

(1998). 

^1 S. Yunoki et al., Phys. Rev. Lett. 80, 845 (1998); E. 
Dagotto ct al, Phys. Rev. B 58, 6414 (1998); A. Morco, 
S. Yunoki and E. Dagotto, Science 283, 2034 (1999) and 
references therein. 

S. Q. Shen and Z. D. Wang, Phys. Rev. B 58, R8877 (1998). 
S. Q. Shen, Phys. Lett. A 235, 403 (1997); S. Q. Shen and 
Z. D. Wang, Phys. Rev. B (to be published). 
A.J. Millis, Nature 392, 438 (1998). 

See A. Auerbach, Interacting Electrons and Quantum Mag- 
netism (Springer-Verlag, Berlin, 1994) 
S. K. Sarkcr, J. Phys.: Condens. Matter 8, L515 (1996). 
S. K. Mishra, R. Pandit, and S. Satpathy, Phys. Rev. B 
56, 2316 (1997). 
2* J. D. Lee and B. I. Min, Phys. Rev. B 55, R14713 (1997); 
L. Sheng and C. S. Ting, Phys. Rev. B 57, 5265 (1998). 

E. L. Nagaev, Phys. Rev.B 58, 2415 (1998). 



7 



